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ABSTRACT 

This  paper  develops  the  concept  of  component  availability  *• 
intrinsically  dependent  upon  the  system  in  which  the  compo^nt  functions. 

The  connection  between  the  repair  aspects  of  the  ion 

availability  is  developed  in  two  ways.  An  integro-differential  equation 
is  derived  that  relates  the  two  concepts.  Alternatively  the  process 

U t«aSd  as  a combinatorial  queue  in  which  each  element  of  the  system  is 
«acS3lhr”s!;  it.  states  of  Availability.  Examples  are  8^ven  in  s^ 
detail  for  small  number  of  elements.  Numerical  solutions  are  given  for 
lu:«r«ir;irpL...  ..th  ti«  v.ryl.«  .»d  .t..dy  .t.t.  «. 

.bwlMd!  S.L  .ppto.cl».  to  . tfd.  oft  .o.ly.l«  for  .v.ll.kUlty 
specification  or  design  are  given. 
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Chapter  U Basis  Concepts  in  Component  Availability 

In  the  study  of  reliability  and  availability  two  major  points  of  view  are 
widely  considered.  On  the  one  hand,  a complete  operational  system  is  con- 
sidered and  various  aspects  of  system  reliability  or  availability  are  addressed. 
Alternatively  one  may  focus  on  the  components  which  make  up  a system  or 
indeed  might  be  employed  in  a variety  of  different  systems.  A widely  con- 
sidered approach  to  these  subjects  is  to  express  system  reliability  or 
availability  in  terms  of  component  reliability  or  availability  by  means  of 
the  laws  of  probability  and  system  logic.  Implicit  in  much  of  this  work  is 
the  concept  of  component  reliability  or  availability  defined  independently  of 
the  system  and  use  for  which  the  component  is  intended.  In  the  case  of 
reliability,  one  can  certainly  give  meaningful,  operational  definition  of 
component  reliability  which  is  completely  independent  of  any  system.  When 
this  is  done,  one  must  still  consider  special  situations  (such  as  stress  and 
shelf -life)  when  applying  specified  component  reliability  to  the  formulation 
of  a system  composed  of  many  parts.  In  the  case  of  availability,  there 
seems  to  be  little  or  no  justification  for  a concept  of  component  availability 
independent  of  the  system  within  which  a component  shall  operate.  The  con- 
cept of  repair  which  distinguishes  availability  from  reliability  carries  with 
it  a number  of  factors  related  to  the  repair  process  which  are  necessarily 
part  of  a system  description. 

This  report  takes  the  point  of  view  that  component  availability 
must  be  defined  within  those  aspects  of  a system  which  relate  to  the  repair 
process.  It  develops  a general  relation  between  component  availability  and 
the  system  repair  process.  It  then  considers  various  repair  process  queues 
in  detail,  using  them  to  obtain  corresponding  component  availability.  System 
availability  is  not  considered  directly  in  this  work,  such  quantities  would  be 
obtained  by  using  the  component  availabilities  as  derived  here  together  with 
system  logic  and  the  laws  of  probability. 

1.  1 Definitions 

In  this  report  a system  is  assumed  to  consist  of  a total  of  N components 
with  specific  components  designated  by  the  index  i=l, . . . , N.  The  study  deals 
with  component  availability  and  does  not  consider  system  logic.  Therefore, 
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component  failures  are  treated  as  events  taking  place  within  the  population 
of  N components.  Several  of  the  quantities  arising  in  the  report  are  defined 
in  this  section,  others  will  be  introduced  as  they  occur. 

Most  of  the  general  material  is  developed  for  a single  "typical"  com- 
ponent. Since  only  that  component  is  involved,  it  is  not  necessary  to 
provide  a subscript  to  quantities  related  to  the  component.  However,  a 
subscript  is  provided  to  clearly  indicate  the  specific  dependence  on  the  com- 
ponent and  to  present  results  in  a form  directly  applicable  to  more  detailed 
investigations  where  more  than  one  component  is  involved. 

Several  of  the  most  basic  quantities  employed  in  this  report  are  now 
defined  and  provided  with  notation: 

ri(t)  = the  probability  that  component  i has  not  failed  up  to  time  t.  This 
is  the  component  reliability  as  a function  of  time.  One  can  consider  a 
"steady  state"  value  for  component  reliability,  defined  as  ^ij^ 

In  most  cases  = 0 however,  for  completeness  it  is  included  among  the 
defined  quantities  arising  in  this  report. 

y.(t)  = the  probability  that  component  i is  operational  at  time  t.  This  is 
the  component  availability  as  a function  of  time.  The  most  basic  aspect  of 
component  availability  is  that  repair  (including  replacement)  is  allowed. 
There  is  no  limit,  in  the  general  situation,  to  the  number  of  repair  events 
that  may  have  occurred  up  to  time  t.  In  fact  the  number  of  such  events  is  a 
random  variable  which,  though  useful  to  study,  falls  outside  the  range  of 
this  report.  The  steady  state  availability  is  well  defined  and  is  most  often 
a positive  quantity  ^iig  yj(t)  = Yj#  though  in  some  special  cases  it  can  be 
zero. 

In  the  initial  work  described  in  this  report,  it  is  useful  to  place  analysis 
in  the  framework  of  exponential  processes.  Such  processes  often  represent 
reasonable  approximations  to  actual  situations  and  at  the  same  time  are 
relatively  simple  to  work  with.  Therefore,  their  use  provides  a meaningful 
starting  point  for  availability  studies.  It  is  interesting  to  note  that  expo- 
nential service  and  exponential  failure  do  not  result  from  an  assumption 
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of  Poisson  arrival  at  the  faile;'  state  from  the  operational  state.  One  requires 
an  infinite  population  of  components,  of  the  same  kind,  in  order  to  make  a 
Poisson  arrival  assumption.  The  assumption  of  exponential  service  (repair) 
time  needs  no  comment.  The  assumption  of  exponential  failure  derives  from 
the  following  general  argument  with  f(x)  a general  failure  density. 

P [a  failure  in  bx\  operational  at  x]  P [operational  at  x] 

= P [a  failure  in  tx  and  operational  at  x]  = f(x)  Ax 

When  f(x)  = Xe"^*, 

Xe"^*Ax 

P[a  failure  in  Ax  I operational  at  xT  — = XAx 

l-[l-e-^] 

The  exponential  repair  and  failure  rates  for  component  i are  denoted  by 
U.  and  X.  respectively. 

The  mathematical  framework  for  the  present  study  is  a queue  which 
represents  repair.  The  total  population  of  N components  consists  of 
operational  components,  failed  components  waiting  for  service,  and 
failed  components  in  service.  Upon  completion  of  service  a component 
returns  to  the  set  of  operational  components.  The  queue  itself  may  be  of 
various  forms  representing  the  repair  features  of  the  system  in  which 
components  under  analysis  are  to  operate.  In  particular  the  queue  will 
involve  specifications  of  service  time  distributions,  number  of  repairmen, 
order  of  selection  of  down  components  for  repair,  and  indeed  the  full  range 
of  possible  variations  to  be  found  in  queueing  models.  This  report  deals 
mostly  with  the  number  of  repairmen  and  specifically  allows  each  component 
to  have  its  own  form  of  repair  time  distribution  (as  specified  by  subscript  i 
in  the  exponential  case  \l^).  In  this  report  the  queue  concept  is  treated  in 
two  distinct  ways. 

‘As  a general  approach  to  relate  the  repair  aspects  of  a system, 
as  reflected  in  the  waiting  time  of  the  queue,  io  component  ' 
availability. 

‘ Detailed  study  of  the  queues  themselves  for  very  small 
(N  s 2,  3)  populations. 

The  waiting  time  distribution  is  basic  to  part  of  this  work,  it  represents 
the  probabilistic  nature  of  the  time  a down  component  must  wait  before 
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service  begins  (in  a broader  context  the  time  to  obtain  a replacement  or 
some  other  part  may  be  included,  with  due  care,  in  the  waiting  time).  The 
waiting  time  is  represented  b'/’  the  distribution  function: 

W.  (s,t)  = P [the  time  component  i waits  in  queue  before  entering 
service  is  < s if  it  joins  the  queue  at  time  t 3.  It  should  be  noted  that  a 
dynamic  version  of  W.  is  defined  where  its  form  is  allowed  to  depend  on  t . 

In  many  cases  W.  will  not  depend  on  t for  one  reason  or  another.  In 
particular,  interest  may  focus  on  the  steady  state  value  |ijn  Wj(8,t)  s W.(8). 

Though  most  of  the  present  work  is  in  terms  of  exponential  service, 
some  use  is  made  of  a general  service  time  distribution  denoted  by: 

■p  (ij  = p [repair  time  for  component  i < tJ,  this  is  not  considered  to 
on  the  time  component  i enters  the  queue  or  enters  service. 

A Fundamental  Relation 

In  this  section  a fundamental  relation  is  developed  between  the  com- 
ponent availability  y . (t)  and  component  waiting  time  behavior  as  represented 
by  W.  (s,t).  That  relation  is  felt  to  provide  a bases  for  a wide  range  of 
investigations  in  the  area  of  component  availability  some  of  which  are 
addressed  in  detail  in  this  report  while  some  others  are  discussed  briefly 
in  the  last  section. 

The  fundamental  relation  is  derived  by  considering  how  component  i can 
be  available  at  time  t+M  from  its  possible  situation  at  time  t.  In  such  a 
development,  it  is  necessary  to  break  the  (possibly  infinite)  chain  of  events 
that  can  contribute  to  availability  at  time  t.  Though  it  is  theoretically 
possible  to  introduce  all  possible  situations  of  failure  and  repair,  it  is 
extremely  cumbersome  to  do  S''.  Therefore,  the  approach  taken  here  is 
to  subsume  all  possible  previous  events  under  the  availability  at  the  previous 
time  values.  The  development  follows. 

y^  (t+At)  is  the  result  of  two  distinct  situations  at  time  t: 

1.  Component  i is  available  at  time  t and  remains  available  over 
the  subsequent  interval  At  (independent  events  by  assumption). 

2.  Component  i is  not  available  at  time  t and  becomes  available 

in  interval  At.  Its  failure  within  At  upon  becoming  available  is  assumed 

of  higher  order  in  At  and  is  omitted  from  the  formulation. 
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If  the  probability  of  situation  (2)  is  denoted  by  H,  then 
y^  (t+At)  = (1-X.At)  yj  (t)  + H 

The  probability  H includes  all  cases  in  which  component  i is  in  the 
service  queue  (unavailable)  at  time  t and  completes  service  by  time 
t+At.  A contribution  is  made  to  H by  every  possible  service  time  for  com- 
ponent i where  that  time  represents  the  service  event  which  prevents  i from 
being  available  at  time  t (and  hence  part  of  situation  (1)  rather  than  H).  For 
component  i to  be  receiving  service  at  time  t,  it  must  have  failed  at  some  previous 
time  T and  entered  the  service  system.  Upon  entry,  it  had  to  wait  a time 
s>  0 (measured  from  t ) and  at  time  t + s its  service  started,  to  be  com- 
pleted between  t and  t+At.  The  possible  cases  for  T can  be  introduced  into 
a sum  of  incremental  terms  which  upon  taking  a limit  process  yields  an 
integral  formulation  for  H. 

Let  A.  = the  probability  that  component  i was  operating  at  time  and 

failed  in  the  next  At  interval. 

A.  = y,  (T.)  X.iT. 

Let  B.  = the  probability  that  component  i spent  the  remaining  time  period 

t+At-T.  in  the  service  system  and  completed  service  it-  ;he  interval  (t,  t+At). 
j m 

B = lim  ( y P [component  waited  till  s,  from  t 

then  entered  service  in  interval  Sj^  + As)  and  from  that  time  on  was  in 

service,  completing  service  in  (t,  t+At)  measured  from  start  of  service  at 

'"i  + ) 

J • n 


Then  B'lg,  I ^ 


The  events  involved  in  the  detailed  formulation  of  Bj  can  be  appreciated 
by  reference  to  Figure  1. 

I Service  I 


t 

meosured  from  zero 


Tj.S 


t-Ti-S 


Figure  I. 


stort 

service 


nAt-Tj-s 


By  definition  of  F.  (t)  and  W.  {s,t)  the  contributions  to  the  probability 
B.  are  given  by  the  following  expression: 

m 

B.  = lim  y CWj{8j^+As,t.)  -W.(Sj^,T.)]CF.(t+At-T..Sj^).F.(t-T.-Sj^):i 
^ k=l 

ni“^ 

where  the  individual  terms  correspond  to  possible  waiting  times  indicated 
by  8^. 

Upon  taking  the  indicated  limits  in  the  expression  for  B^,  one  obtains  the 
Stieltjes'  integral: 

tk  t-T. 

B.  = J CFjJ(t+At-T^-8)  .Fj(t-T.-s)  dWj(8,  T^). 

Introduction  of  this  expression  into  the  equation  for  H and  form\ilation  of 
the  limits  in  H yields; 

H = I yj(T)  J[Fj  (t+At-T-s)  .Fj  (t-T-8)T  dWj  (s,  t)  dr. 
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Returning  to  the  formulation  of  the  equation  for  availability  of 
component  i and  dividing  through  by  At  one  obtains: 


y.  (t+At)-y.  (t) 


t-T 


“ST 


= -X.  y.(t)  + r y.  (,)  X.  J {F.  (t+6t.8)  - F.(t-T-s)l^„  ,^ 


In  the  limit  as  At-«,  under  the  assumption  that  the  service  time 
distribution  possesses  a density  function  L (t),  the  availability  relation 
becomes: 


dyAt)  • . 

y (t)  + I y (T)X  1 f (t-T.s)  dW  (s,T)  dT  (1) 

Equation  (1)  is  the  fundamental  relation  that  was  desired.  It  relates 
component  availability  y^it)  to  the  waiting  time  distribution  W^(s,t)  for  the 
same  component.  The  waiting  time  represents  the  system  in  which  the 
component  is  to  operate.  Features  of  the  other  components  and  of  the 
system  repair  and  service  policies  are  all  contained  within  the  W.(s,t) 
expre  ssion. 

Equation  (1)  utilizes  the  Stieltjes'  integral  employing  waiting  time 
distribution  W.(s.t)  because  that  distribution  often  has  a finite  jump  at  the 
8=0  value  which  presents  the  general  use  of  a density  function  representation. 

In  this  section,  two  extreme  cases  for  equation  (1)  will  be  evaluated, 
an  additional  case  will  be  illustrated  in  section  2.  1.  3.  Potential  applications 
and  investigations  involving  equation  (1)  are  discussed  in  Chapter  6. 

Case  1 no  repairmen.  This  is  one  extrem.e  case  in  the  concept  of 
component  availability.  It  means  that  a component  is  only  available  until 
it  breaks  down.  Thus  availability  is  the  same  as  reliability  for  this  case. 

The  waiting  time  for  repair  is  infinite  so  that  W^(s,  "tIsO  for  all  finite  values 
of  s.  Thus  the  Stieltjes  integral  over  W.(s,  t)  is  zero  and  equation  (1)  becomes: 

dy  . 

■ar  ""  "^i  ^i 
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which  is  the  well  known  equation  for  component  reliability  with  solution: 

y.  (t)  = r.  (t)  = 

Case  2 repairman  for  each  component.  This  is  the  extreme  case 
opposite  to  Case  1.  In  this  case  each  component  has  a repairman  available 
to  start  repair  whenever  that  component  breaks  down.  Thus  waiting  time  is 
always  zero  so  that  W.  (s.t)  =1  for  s > 0 . The  Stieltjes  integral  con- 
tributes a single  term  at  s=0  and  gives  zero  values  elsewhere. 

It  is  convenient  to  specialize  this  case  further  by  assuming  exponential 
service  time  distribution  with  service  rate  u..  These  observations  yield 
the  following  form  for  equation  (1): 


dy  .(t) 


^ . -X.y  . (t)  + U;  X.  .-“i‘  [ Yi  (T)  eW"  dT, 


or  alternatively; 


f d y . (t)  1 


UiT 

(t)  e dT. 


Differentiaticnwith  respect  to  t (using  Leibniz's  rule  on  the  right 
hand  side)  yields: 


d y.  (t) 


dt 


d y.  (t) 

+ ST' 


= 0. 


+ k. 


with  general  solution. 

ttSt 

It  remains  to  evaluate  the  constants  c and  k.  At  t=0  it  is  assumed  that 
y.  (0)  = 1.  To  obtain  another  relation  one  observes  that  near  t=0,  before  the 
first  failure  occurs  availability  is  the  same  as  reliabUity.  Thus  it  is  assumed 
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that  these  quantities  have  the  same  slope  (first  derivative)  at  t=0.  Since 


dr.(t) 


drj(t)  I ^ I r>r  r=-\..  One  then  obtains  k=  v— . 

This  leads  to  the  well  known  result  (reference  l,page  337)  for 
availability  when  repair  is  undertaken  at  once  with  exponential  service 
time  distribution: 


= -X. 


. e’^i*  ard  = ce'^^i^^i^  ^ the  assumption  yields: 


dyjt) 


yi(t)  , 

with  steady  state  value 


X.  + Ui 
1 1 


.(X  + Uj)  t 

e 1 I + 


Yi  - t-« 


[-Tu: 


Though  this  is  a well  known  result,  it  is  derived  here  in  a different  way 
than  elsewhere  in  the  literature.  It  follows  from  and  illustrates  equation  (1). 

In  this  context  it  is  very  clear  how  the  availability  and  the  system,  represented 
by  waiting  time  for  service,  are  related. 

The  remainder  of  this  report  deals  primarUy  with  studying  the  detaUs 
of  the  system  service  procedure  which  determines  waiting  time  distribution 
W.(s,t).  Such  considerations  also  yield  availability  directly.  In  Chapter  6 
some  further  consideration  is  given  to  the  interplay  between  availability 
and  waiting  time  as  expressed  by  equation  (1). 

1.  3 Combinational  Queues 

This  report  deals  with  component  availability  within  the  operational 
context  of  a system  comprised  of  a number  of  components.  The  system 
has  some  facility  to  repair  (or  replace)  components  which  provides  the, 
actual  level  of  component  availability.  A relation  between  the  system 
service  capability,  as  expressed  by  component  waiting  time,  and  the 
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component  availability  has  been  described  above  and  given  mathematical 
form  in  equation  (1).  More  detailed  study  of  the  waiting  time  and  the 
availability  may  be  carried  out  by  means  of  a queue  model. 

It  is  assumed  that  a system  consists  of  mcomponents  and  that  each 
component  is  characterized  by  its  failure  rate  and  its  repair  rate  u.. 

More  general  probability  forms  can  be  used  but  are  difficult  to  work  with 
(short  of  simulation)  and  this  report  limits  itself  to  exponential  processes. 
The  system  is  provided  with  a service  capability,  specified  in  each  case 
under  study,  which  can  be  generally  referred  to  as  the  service  system,  or 
service  queue.  When  a component  fails,  it  enters  the  service  queue,  waits 
if  necessary  for  service,  receives  service,  and  upon  completion  of  service 
returns  to  operation  in  the  system.  The  history  of  failure  and  repair  for  a 
component  represents  the  component  availability.  In  this  process  of  com- 
ponent failure,  repair,  and  operation  the  system  logic  plays  no  role.  From 
a broader  point  of  view  the  interaction  of  components  in  system  operation 
may  effect  individual  properties  such  as  failure  rate  X but  any  such  con- 
siderations are  assumed  to  be  accounted  for  by  parameter  values  so  far 
as  component  availability  is  concerned.  System  availability  must  consider 
the  combination  of  component  availabilities  following  the  system  logic.  The 
system  itself  and  the  set  of  all  components  effect  individual  component 
availability  only  through  the  service  queue. 

Thus  one  deals  with  a queue  system  having  a finite,  fixed  total  population 
of  m district  elements,  each  of  which  has  its  own  parameter  values  (proba- 
bilistic behavior)  in  general.  There  are  three  ' states"  possible  for  any 
member  of  the  population:  operational  (available),  waiting  for  service,  or 
in  service  (the  last  two  correspond  to  the  unavailable  condition). 

Such  a model  can  be  thought  of  as  a cyclic  queue,  a topic  that  has 
received  attention  in  the  literature  (for  example  in  reference  2).  However, 
the  interest  in  this  report  is  to  a special  kind  of  detail  and  result  not  con- 
sidered in  the  standard  approaches  to  cyclic  queues.  For  the  kind  of  results 
desired  here, it  is  necessary  to  follow  each  member  of  the  population  through 
the  queue  process.  For  m of  any  size,this  becomes  extremely  involved  to 
the  point  of  practical  infeasibili^.  However,  for  m=2  and  3,  it  is  completely 


JL 
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feasible  and  is  done  in  the  following  chapters. 

Because  of  the  detailed  nature  of  the  queues  treated  here,  it  seems 
appropriate  to  give  them  the  name  ' combinatorial  queues''.  This  implies 
full  attention  to  the  individual  character  of  each  element  of  a (small)  finite 
population  as  it  interacts  with  the  other  members  through  the  queue  model. 
The  notation  used  for  these  models  is  form  elements  and  r repairmen 

(channels).  In  each  case  first  come  first  served  (FIFO)  service  discipline 
is  assumed  and  none  of  the  possible  variations  of- queue  behavior  are  con- 
sidered. Many  generalizations  suggest  themselves  for  future  study.  An 
important  special  case  is  when  all  the  components  (population  elements) 
are  the  same,  i.e.  X.  ad  do  not  depend  on  i.  For  these  cases  the  special 

notation  Q®  may  be  used, 
m,  r 

The  cases  for  m=2  andm=3,  though  small  size  from  a general  point  of 
view,  are  studied  in  detail  for  two  reasons: 

• These  cases  are  small  enough  to  be  given  through  investigation  indicatiiig 
a variety  of  features  in  the  nature  of  combinatorial  queues  as  they  relate  to 
component  availability. 

• Components  may  be  relatively  complex  and  significant  parts  of  a system 
so  that  actual  systems  may  be  encountered  which  are  usefully  considered  as 
having  only  a few  components.  For  such  systems  detailed  studies  as  illus- 
trated for  the  m*2  or  3 cases  are  appropriate.  When  this  is  the  case,  it  is 
felt  that  the  combinatorial  queue  approach  is  likely  to  be  the  most  powerful 
treatment  of  system  availability.  By  such  methods,  a complete  analysis 
(without  simulation)  can  be  produced.  The  kind  of  results  and  their  use  is 
indicated  by  what  is  done  for  the  cases  in  the  next  two  chapters  ahd  also  in 
Chapter  5. 


Chapter  2.  Two  Component  Queues. 

This  chapter  deals  with  queues  of  the  form  Q2  The  cases  r«0,  no 
repairmen  and  r=2,  full  number  of  repairmen  have  been  given  in  Section 
1.  2 where  componen^  availability  is  obtained  from  equation  (1)  for  these 
cases  in  general.  The  only  case  to  be  considered  here  is  whenrsl,  one 
repairman  to  service  the  components  as  they  fail.  Though  one  could  form- 
ulate the  Q queue  analysis  in  terms  of  Markov  chains  and  the  states 
V 
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of  the  system,  the  approach  of  this  report  is  to  use  differential 
equations  as  the  mathematical  representation  of  the  queue  system. 

2.  1 One  Repairman 

The  notation  of  this  section  is  as  follows: 
y.  (t)  = availability  of  component  i at  time  t,  i=l,  2 . 

= steady  state  availability  of  component  i. 

P..  (t)  = probability  that  at  time  t component  j is  in  service  and 
component  i is  waiting  for  service.  Only  certain  values  of  i and  j lead  to 
possible  forms  of  namely  Pqq,  Pqj,  Pq2*  ^12»  ^21* 
zero  indicates  no  component  in  that  situation.  For  example  P^p  means 
both  components  are  operational.  When  steady  state  values  are  intended, 
the  same  notation  is  used  without  the  variable  t. 

In  this  notation,  availability  is  the  condition  of  not  being  in  the  service 
system  so  that: 

y,  (t)  . Poo  (t)  + Poj  (1), 
y J (t)  = Poo  (t)  + (t). 

Waiting  time  can  also  be  expressed  in  terms  of  the  queue  probac. ilities 

P..  (t)  or  p ,,  in  the  steady  state.  For  simplicity  consider  the  steady  state 

waiting  time  distribution  W.  (s)  = W^  (0)  + Pr  [ a completion  in  < s I arrival 

finds  one  in  servicel.P  , where  P =Pr  [on  arrival  one  is  in  service].  Thus 

^ ft  ft 

one  must  consider  the  a posteriori  probabilities. 

Pr  [component  i completes  service  ^s  | another  customer  finds  i in 
service]  = l-e”^i®  by  the  assumption  of  exponential  service  time. 

Let  q^Q  (1)  = Pr  [component  1 arrives  and  finds  none  in  the  service  system]. 
qQQ  (2)  a R: [component  2 arrives  and  finds  none  in  the  service  system] 
qQ2  = Pr  [on  arrival  component  1 finds  2 in  service] 
q . = Pr  [on  arrival  component  2 finds  1 in  service] 


Then: 


Wj(s)  a q^Q  (1)  + ( Ue  -^^2") 

W2  (.)  * qoo  (2)  + (Ue  -^1*)  q^j 


These  values,  upon  some  simplification  yield  the  following  expressions 
for  waiting  time  distributions: 


The  same  logical  cases  contribute  to  the  general  time  dependent  waiting 
time  distributions  so  that  one  has  the  following  expr'Sssions; 


2.  1.  1 Differential  Equations 

In  the  subsections  of  this  section,  the  differential  equations 
for  the  five  service  system  probabilities  Pjj  (t)  are  presented.  The  first 
subsection  gives  the  general  case  and  the  next  subsection  gives  the  case 
when  both  components  are  the  same.  Logical  cases  contributing  to  the 
probabilities  are  combined  to  form  the  equations,  in  most  cases  terms  in 
(At)*  or  higher  are  never  shown  since  they  drop  out  upon  forming  differential 
limits.  It  is  felt  that  the  various  cases  are  clearly  indicated  within  the 
equations  themselves  and  further  discussion  is  therefore  omitted. 

2.  1.  1,  1 General  Case 

The  initial  formulation  for  the  combinatorial  queue  | gives  the 
following  system  of  equations: 

Pqq  (t+At)  = (1  - XjAt)  (1  - >2At)  PQQ(t)  + (1  - >2^^)  UjAt 

+ (1  - >jAt)  U2At  Pq2  (t) 

Pqj  ( t+At)  = (1  - (1  - UjAt)  Pqj  (t) 

+ >lAt  Pqq  (t)  + U2At  Pj2  (t) 

Pq2  (t+At)  = (1  - XjAt)  (1  - Pq2  (t) 

+ )2At  pQQ  (t)  + UjAt  P21  (t) 

Pj2  (t+At)  = At  (1  - U2At)  Pq2  (t)  + (1  - U2^)  ^13  (t) 

P^j  (t+At)  = X^At  (1  - UjAt)  Pqj  (t)  + (1  - UjAt)  P2,  (t) 

Upon  renewal  of  some  higher  order  terms  in  At  the  system  becomes: 

Pqq  (t+At)  = (1  - >jAt.X2At)  Pqq  (t)  + UjAt  Pqj  (t)  + U2At  Pq2  (t) 

Pqj  (t+At)  = (1  - - UjAt)  pQj  (t)  + X,  At  Pqj  (t)  + a2At  Pj2  (t) 

Pp2  (t+At)  = (1  - XjAt  - U2At)  Pq2  (t)  + ^21^^^ 

Pj2  (t+At)  = XjAt  Pq2  (t)  + (1  - U2At)  Pj2(t) 

P2J  (t+At)  a X2At  PQj(t)  + (1  - UjAt)  p2i(t) 
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Forming  differential  quotients  and  taking  the  limit  as  &t-*0  results  in 
the  following  system  of  difference  differential  equations.  Since  the 
difference  is  accounted  for  completely  by  explicit  expression  of  each 
function,  the  system  is  in  fact  a system  of  differential  equations. 


dPoi(t) 


= -(XjtXj)  Poo'*'+“l  Pol'**  * “2  Poz'*' 

= -'»2«'l>  Poi**''*‘lPoO<*>'X“2Pl2<*> 


. Po2<«*'2Poo“>  + “l  Pzi**' 


dP^,(t) 


" "^1  ^21**^ 


One  can  in  principle  solve  this  linear  system  by  transform 
or  other  "elementary"  methods.  Such  solutions  become  extremely  com- 
plicated. Therefore,  it  is  felt  that  numerical  solutions,  for  selected 
parameter  values  should  be  used  to  obtain  desired  solutions  and  to 
appreciate  the  nature  of  such  solutions,  some  such  results  are  given  in 
the  Appendix.  Another  approach  is  to  study  the  steady  state  as  done  in  a 
subsequent  section.  Still  another  approach  is  to  consider  the  similar 
components  case  Q®  , which  is  done  in  the  next  subsection.  Before  turning 
to  that  topic,  the  "transformed"system  is  expressed  below  as  a matter  of 
general  reference.  To  formulate  that  system,  the  notation  «sed 


for  the  (Laplace)  transform  of  Pij(t)  and  assumed  initial  values  are  P.j(0)=0 
except  for  i=j=0  and  ^00(0)=^*  results  in  (the  dependence  on  s is  not 

explicitly  shown): 
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(a+>  ^+>3)  Pqo 


■'-‘iPoi 


•^2^02 


-^POO  + POI 


■^2  Poo 


■^'zPiz 


+(s+)j+U2)  P02 

->lP02  <s+U2)Pi2 


-^iPzi 


+('3+Ui)P21 


2.  1.  1.  2 All  Components  Are  the  Same 


In  this  case  > 1=^2”^  ^01^*^  ~ ^02^^^ 

Pl2(t)  = ^21^*^'  typical  values  PQjlt),  ^12^*^^  transforms 

P,  2(3)  will  be  used.  There  are  now  only  three  quantities  to  be 
determined  rather  than  five.  The  transformed  system  is: 


(s+2^)  Pqq  -2u?fQl 

->ffoo  9oi  “P12 

-’Poi  ■'■••♦“•I’lZ 


From  which  one  obtains: 


TTTT 


« + (2u  * X)  . 

Poo'”*  - . (.-r,)  (t-rj) 


where 


-(3X  + 2u)  4-  4)m 


rj= ^ 


^2  = 


_ -(3) 


+ 2U)  -J 


+ 4>u 


Standard  algebraic  reduction  and  partial  fraction  expansion  yields  the 
solution  in  terms  of  the  quantities: 


1 

a = — 2" 

2X  + 2>u  + U 


b = 


‘Ju  + 4TB 


c = 


The  results  are; 


Poo(t)=U^a+(ri^+r(2u+>)+U^)b  e^l*+  lr2^+r2(2u+X)+u^)c  e 

P j(t)  = U)a  + X(rj-Hi)  be'l*+  X (r2-Hi)  c 


P02(t)  = Poj(t) 

Pj2(t)  = >^a  + >^be'’l*+x2c 


P2i(t)  = ^12 

The  values  r j and  r2  are  negative  for  all  positive  values  of  X and  u 
(which  are  never  negative,  being  rates  of  events).  Thus  limit  values  exist 
for  all  the  P..(t)  quantities  for  all  X.ji values.  In  infinite  population  queues 
this  is  often  not  the  case  but  here  all  the  probabilities  are  defined  for  any 
time. 

Direct  calculation  establishes  that  the  sum 


Pqq  (t)  + 2 PqjW  + 2 * I as  it  should  be. 

The  steady  state  values  when  the  components  are  the  same  rnay  be 
found  from  the  above  expressions  to  be: 


s lim  Pftft(t)  = U a = 


00 


1 

1+2  {^)+  2(i)2 


PqI  " PqZ  ■ ^00^^^ 


Pl2  “ P2I  “ ^ U ^ PQO 

t-«6 

This  result  will  also  be  obtained  as  a special  case  of  the  general  steady 
state  results  in  Section  2.  1.  2.  2. 

2,  1, 2 Steady  State 

s 

In  this  section  the  steady  state  values  for  j and  j are  obtained. 

This  is  done  by  supposing  that  the  probabilities  have  become  constant  over 
time  so  that  their  derivatives  are  identically  zero.  A system  of  algebraic 
equations  is  obtained  for  the  steady  state  probabilities.  Lower  case  letters  are 
used  so  that  p..  is  the  steady  state  version  of 

2.  1.2.  1 The  General  Case 

The  equations  presented  in  Section  2.  1.  1.  1 yield  the  following  steady 
state  system: 


.(1 

l+>2) 

Poo  + ^^1 

Poi 

+ 

^^2 

P02  " 

-(X 

2‘*‘^1^ 

Poi  + h 

Poo 

+ 

^2 

Pl2  = 

-0 

l+^2> 

P02  + ^2 

Poo 

+ 

^1 

P2I  = 

- a 

2 P12 

+ \ P02 

= 0 

-^1  P21 

'2  Poi 

= 0 

Solution  is  (by  successive  substitution) 

P02  “ ^2  >,  u,  + a 


pi  + >2  ^1 1 

[>j»il  + >2U2  + Ui  U2j  ^00 


) j + X2  + U2 


1 Ui  + >2  U2  + UjUz 


Poi=  h 


Pzi  = 


>1^2 


U 


1 


+ >2  + U2 

>1  Ul  + >2  ^2  + Ui  U2 


Poo 


since  the  sum  of  the  probabilities  must  be  unity,  we  get  by  forming 

1 


that  sum; 


S = 1 + A. 


Poo  = 5 


where 


’X1+X2+U1 

'1^2 

> i+^2'’’^l 

^2 

S — 

+ Xj 

A 

^2 

A 

■ ■ 

^1^2 

U, 


Xi+>2+U2 

—J. 


, where  A=X JUJ+X2U2+WJU2 


The  expression  for  S and  indeed  for  any  of  the  above  values  does  not 
appear  to  simplify  for  general  values  of  Xj,  X2f  Uj>  3-nd  (i2.  To  study  them, 
one  must  employ  numerical  evaluation  one  form  of  which  is  discussed  in 
Chapter  5.  Other  numerical  results  are  presented  in  the  Appendix. 

2.  1.2.2  Components  The  Same 

When  Xj  = X2  = > and  = U2  = U are  introduced  into  the  general  steady 
state  solution  given  in  the  previous  section  one  obtains  directly: 

S = 1 + 2 (^)  + 2(^)2 


and 


Poo 


l+2(^)  + 2{i)  ^ 


i\)i  " P02  " y Poo 

,X.2 

P12  " P21  " 'w^  Poo 

in  agreement  with  the  values  previously  obtained  in  section  2.  1.  1.  2. 


2.  1.  3 An  Illustration  of  the  Fundamental  Relation 

As  described  above,  the  j results  can  be  used  to  formulate 
y.(t)  and  w.(s,t)  for  this  case.  It  is  interesting  to  introduce  these  values 
for  one  component,  say  i=l,  into  the  fundamental  relation,  equation  (1)  of 
Section  1.2  and  observe  that  the  relation  does  indeed  yield  an  identity  under 
such  substitution.  This  gives  a kind  of  closure  to  the  two  distinct  approaches 
developed  in  this  report  and  indicates  the  complimentary  role  of  each  approach 
in  the  study  of  component  availability. 

Set  i=l  in  equation  (1)  of  Section  1.2  and  assume  exponential  service 
time  distribution  to  obtain: 


dyj(t) 


>iyi(t)  + Jyi('r)>i  d Wj(s,  t)  dT 


-w 

Let 


j,  [*Vi  dWj(s,T),  then  using  material  introduced  previously 


for  the  waiting  time  distribution  one  obtains: 


J = U,  w,(0,  t)  + fui  e ^ ^2^02^’’^^ 

A ^ r\  .lu."  r'  I . ^ 


."UZ®  ds 


where  d Wj(s,  t) 


ds 


Po0<’'+‘'02*"’ 


-W28. 


w,(0,  T)  = 


P00'’>+’’02<'''  ■ 

Simplification  and  evaluation  of  the  integral  in  J leads  to: 
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From  previous  work  yi('^)=PoO^'^^'^^02^'^^  ®°  present  version  of 
equation  (1)  becomes; 


= -'itPooC'-^Poz*''*! 


>^02^*^^  +W2Po2^^^  e-^'2(t-^)jdT 


Alternatively: 


<i[P„n(’')+Po2l''>l 

dt  -‘lfP00<'''+P02<’>l 


^r\ 


+ Ul>j  e-  1 / Poo('r)  e 


^1*^  dT4-  P„,(t)  e^2'^dT 


ai-U2 


_ }V^  e'^lt^V  (t)  e^^l"^  d^}  . (a) 

U1-U2  J 02  J 


Turning  to  the  queue  equations  for  Q2  ^ shows  that: 


d[Pnn<'^)+Pn,(’^)] 


00'  ' 02' 
dt 


' (^1  + ^2^  ^00^^^  '*’  ^1^01^*^  ■*’  1^2^02^^^ 
- (>i  + \X-^  ^02^*^  ’’’  ^2^00^^^  '*’  ^1^21^^^ 

->l[Poo(t)  + Po2<*)l  + + 


Therefore  the  present  form  of  equation  (1)  will  be  established  as  an  identity 
if  the  bracket  expression  in  equation  (a)  is  shown  to  be  Wj[Pq j(t)+P2j(t)]. 


\ 
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For  simplicity  the  Hj  multiplier  may  be  divided  out  at  once  since  it  is 
present  in  every  term  of  the  bracket. 

From  the  differential  equations,  one  obtains: 


d[Poi(t)  + P2,(t)l 


where  d 


, -Ui[Po,(t)  + PjiWi  t ’.Poo*"  + “2^12*" 


so  that  ^ 1 ® 


>tr%U2t 


^02^^^  dT  which  may  be  introduced  into 


the  differential  equation  above  to  give 


^ o 

which  upon  integration  yields: 

Poi(t)  + P2^(t)  = >1  Poo(w)  dw 


+ U2> 


1 


jw  p (T)dT  dw 


The  first  term  on  the  right  is  equal  to  the  first  term  in  the  bracket  expression 
in  (a)  with  Uj  divided  out.  Integration  of  the  second  term  by  parts  yields  the 
remaining  two  terms  of  the  bracket  expression  as  follows: 

Let  u=r  e^Z"^  Pg2(T)dT  and  dv  = e^^^l'^Z^^dw  then  the  technique  of 

X) 

integration  by  parts  yields; 

il2>j  e'^l‘  I*  e^Z"^  Pq2('^)  d^dw  = 


n. 


r 


I ! 


n 


I i 


it 


I ! 


I 


?. 

y ! 


liir 
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e'^1^ 


(Ul-U2)t  UoT 


/ 


>“2  PozW'^-'- 


Y Poz*"'  '*■' 

< “l-“2 


u,-u- 


PojindT 


‘1-2 


U2),e-“l'  re“l'' 


PojlTjdT 


as  required. 

This  illustrates  that  introduction  of  the  values  w.(s,t)  and  y^(t)  obtained 
from  the  solution  of  Q,  , satisfy  the  fundamental  relation,  equation  (1), 


as  in  identity. 

2.  2 Summary  Results. 

Because  of  the  relatively  complex  form  of  the  solutions  even  in  the  two 
component  case,  general  results  require  numerical  calculation.  Some  such 
considerations  are  presented  in  Chapter  5 and  the  Appendix.  In  this  section 
the  case  in  which  both  components  are  the  same  is  summarized  in  terms  of 
the  algebraic  solution  forms. 


s s 

For  two  components,  the  only  cases  possible  are  Q2,  0,  Q2 


‘2.  2 


so  long  as  a FIFO  ervice  discipline  is  assumed.  This  report  does 


not  consider  generalizations  of  the  combinatorial  queue  problem.  Q2  q has 


the  (reliability  case)  solution 
' t 

y(t)  = e“  for  component  availability. 


The  steady  state  value  ^ y(*)  = ® ca8e.Q2  jhas  solution 

y(t)=(u^+aX)a+[rj  + r j(2u+X)+u^+X(rj+u)]b  e^^^+Cr^  +r  2(2u+X)+U^+>  (r2+U)]c  e’‘2* 


whe  re  ^ (3>  ■t2u) 


+ 7x^  + 4Xu  . _ (3X  + 2u)-V  X^  + 4Xu 

— , To-  - 9 


a= 


1 


2X^  + 2>u+ 


2 

, b= 


c= 


2 

-1 


4Xu 


4X^ 


1 

f , 

M 


r 


\^p 


q|^ 


Q2,0 


Figure  2 

Chapter  3 Three  Component  Queues 

This  chapter  deals  with  a system  having  three  components.  The  FIFO 
assumption  is  made  as  is  the  assumption  of  exponential  service  and 
exponential  breakdown  distributions.  Special  cases  for  steady  state  and  all 
components  the  same  are  considered.  There  are  a total  of  four  possible 
kinds  of  repair  queue:  q,  Q3  i»  ^3  2’  ^3,  3'  first  and  last 

of  these  have  already  been  treated  in  the  general  results  of  Section  1,2  . 

This  chapter  treats  j and  2*  Results  become  complicated  and  can 
not  readily  be  carried  as  far  as  the  two  component  cases.  It  is  fait  that 
for  more  components  special  studies  are  called  for  and  simulation  of  cases 
of  interest  might  prove  as  (or  more)  effective  as  the  differential  equation 
approach.  However,  thp  material  in  this  chapter  should  serve  as  a guide  to 
the  more  involved  cases.  It  is,  of  course,  useful  in  its  own  right  when  a 
system  consists  (in  some  operational  configuration)  of  three  components 
(likely  to  be  individually  complex). 

3.  1 One  Repairman 

In  this  case  the  notation  Pjjj^(t)  is  used  for  the  probability  that  at  time  t, 
component  k is  in  service,  j is  next  in  line  waiting  for  service,  and  i is 
second  in  line  waiting  under  the  FIFO  discipline.  The  value  0 means  no 
component  occupies  the  indicated  position.  Possible  service  queue 
probabilities  are; 


¥ 


i oiniii  i»w.  m iwi^ii  I 
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Pqoo  meaning  none  in  the  service  system, 

Poor  P002’  P003’  system  and  in  service, 

P02l'  P031'  P012'P032*P013'  P023'  one  is  waiting  and  one 

in  service, 

Pl23’  P2I3’  Pl32’  P3I2'  ^231’  P32!/  two  a--' one 

in  service. 

Availability  is  given  by  appropriate  combinations  of  the  service  system 
probabilities,  for  example: 

yj(t)  = P*000^^^  ^002^^^  ^003^^^  ^032^^^  ^ ^023^*^ 

(i.e,  component  1 is  not  in  the  service  system).  Waiting  time  can  also  be 

considered  as  illustrated  in  Chapter  2. 

3.  1.  1 Differential  Equations 

In  formulating  the  differential  equations,  higher  order  terms  in  tt  are 
either  not  included  or  are  removed  as  limits  are  taken.  The  equations  are 
formed  by  considering  the  logical  cases  that  contribute  to  the  probability 
values  under  the  exponential  distribution  assumptions; 

3.  1.  1.  1 The  General  Case 

The  system  of  equations  for  queue  probabilities  at  t+At  in  terms  of 
values  at  t follows: 

1 -(>  j+>2+>  3)  At]  P*000^^^  ^1^^  ^00  ^002^*^'''^3  ^00 3^^^^ 

Pooi<t+At)=ri-UjAt]ri.(>2+>3)At]  Pooi(t)+>iAt  Pooo(*)+^^2^^  Poi2<*>+^"3^*^013<*) 

Poo3(t+  At)=ri  -U3  At][  1-(X  j+>  2}  At]  ^003^*^'*’^  ^031  ^*^■’’^2 ^*^032^^^ 

P02i(‘+At)=(l-Uj  At){I-)3At)Po2i(t)+>2^  ^001^*^  ^ ‘^3*‘  ^213^*^ 

P03i(t+At)=(l-Ui  At)(l->2At)Po3i{t)  + > PqOI^*^  ■*■  ^^2^^  ^312^*^ 
Poi2(t+At)=(l-U2At)(l->3At)Poi2(t)+>iAt  Pqoz^^]  +^3^*  ^123^*^ 
P032(t+At)=(l-U2At)(l.)  jAt)Po32(t)+>3At  Poo2(t)  + UiAtP32i(t) 
Poi3(t+At)=(l-U3At)(l->2At)  PQ^3(t)+>j&t  Poo3(*]+  *^2^*  ^132<^^ 
P023(^+At)=(l-U3At)(l->jAt)  Po23(*)+^2^t  l^l^*  ^231^‘^ 

Pj23{t+At)=(l-U3At)  Pj23(t)  + ^023<*^ 

P2j3(t+At)=(l-U3At)  P213  (t)  + >2At 


Pj32  (t+At)=(l-U2At)  Pi32(t)  + >j&t  PQ32(t) 

P3I2  (t+frt)=(l-U2frt)  P3i2(t)  + ^ 

P23J  (t+At)=(l-Uj&t)  P23i(t)  + >2^* 

P32J  (t+M)=(l-Ujfrt)  P32i(t)  + >3^  Po2i(t) 

Upon  division  by  At  and  passing  to  the  limit  as  &t-*0,  one  obtains  the 
following  system  of  differential  equations: 


**^00 1 


dPo2i(t) 


dPo3,(t) 


dPoi2(t) 


■iPosz'** 


dPoj3(t) 


dP,23(t) 

3t 

dP2i3(t) 

3t 

dPi32(t) 


-(>  ^+>2+^3)  '^1^001^*^  ^2  ^002^*^  ^^3  ^003^*^ 

-(Uj+>2+^3)  ^001^^^  ^ ■*■  1^7  ^017^*)  ^ni 


r 000'-'  ^2  "012'-'  ■ ^*3  * 013' 


-(U2+^1+>3)  ^002^^^  '*'  ^2^000^*^  ^1  ^021^*^  '*’ 


3 "023' 


-(U3+>l+^2^  PQ03(t)  + ^3^000^*^  '*’  ‘^1^031^*^  '*'  ^2  ^032^*^ 


•(aj+>3)  Po21^*^  ^2  ^001^*^  ^3  ^213^*^ 


-(U2+^2^  ^031^^^  ^ ^ 


3 "001'-'  " ‘ 312' 


■(a2+'3)  ^012^*^  ■’■  ^1^002^*^  ^3  ^123(t) 


-{U2+>i)  ^032^*^  ■*'  '3  ^002^*^  ■’■  ^1  ^321^^^ 


-(U3+)2)  ^013^*^  ■*■  ^1  ^003^^^  ■*■  ^2  ^132^*^ 


.(U3+>2)  PQ23<t)  + ^2  ^003^*^  * ^1  ^231^^^ 


■*^3  ^123^*^  ^1  ^023^^^ 


■^3  ^213^*^  '*'  ^2  ^013^*^ 


■^2  ^132^^^  * ^1  ^032^^^^ 


dpjljd) 


dP23,(t) 


“ ■'^2  ^3 


— -dt “ ‘^1  ^^231^*^  ■’■  ^2  ^031^*^ 


dP32i(t) 


“ "^1  ^321^*^  '’’  ^3  ^021^^^ 


This  system  can  only  be  treated  by  numerical  solution  methods.  Some 
results  are  presented  in  the  Appendix,  The  steady  state  can  be  considered 
along  with  some  special  cases.  This  is  done  in  the  present  chapter. 

3.  1.  1.  2 All  Components  The  Same 

When  the  three  components  are  the  same  yj=U2=U2=U  and  ^ 

The  sixteen  probabilities  reduce  to  four,  which  are  taken  to  be  (as  typical 
cases)  Pooo(t).  Pqoi^*)'  ^012^*^  Pl23<*)-  differential  equations 

become; 


dPnnn(t) 


= -3>Pooo(t)  + 3uPooi(t) 


dPn,,(t) 


= .(U+2X)  Pooi(t)  + > Pooo(t)  + 2u  Poi2(t) 

= -(U+»  Poi2<*)  + ^ ^001^^)  + ^123<^^ 


dPi,-(t) 

. -u  Pj„(i)  + > P„,2W 

Though  considerably  simpler  than  the  full  system  when  the  components 
differ,  this  set  of  equations  does  not  seem  to  yield  a simple  closed  form 
solution  with  reasonable  effort.  One  approach  to  solution  is  to  take  the 
Laplace  transforms  using  initial  conditions  other 

probability  values  zero  at  t=0.  This  results  in  the  following  set  of  equations 
in  the  transform  variable  s. 

• ^ooo<*>  - '■ooo"»  = Pooo'’>  ^ 

• Pool'*’  - Pooi'°>  = Pool'*’  + ’ Pooo'*’  + ^ Poi2<*’ 
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* ^012^®^  " ^012^°^  ~ -(U+X)  ^012^*^  ■*■  ^ ^001^®^  ''"  ^ ^123^®^ 

® ^123^®^  " ^123^°^  " ^123^®^  ■’■  ^ ^012^®^ 

which  become,  upon  substitution  of  initial  values: 

(s+3>)  Pqoo^®)  = 1 + 3u  PQOI^®^ 

(s+u+2>)  PqOI^®^  " ^ ^000^®^  ^012^®^ 

(s+y+X)  P ^ ^001^®^  + U P j23^®^ 

(s+u)  Pi23<®)  = ^ P012(®^ 

It  should  be  kept  in  mind  that  there  are  in  fact  6 values  like  Pj23»  ^ 

like  Pqj2  ^ ^001'  being  representative  values. 


The  above  system  yields; 
^123^®^  "81^  ^012^®^' 


^s  ^001^®^ 


Poi,(®)  = - ^■■- Pnni(®) 


(s^•^)  +ls 


^001^®^ 


>[(s+u)^  + Xs] 

= 

(s+u)  + 3>8  (s+u)  + 2>  s 


^000^®^ 


So  that: 


P,,,(8)  = 


123''"'h7sr  ^000 


Pnnn(®) 


^012^®^  ■ h/s)^^^  ^000*®^ 


7s  ,.s  _ >r{8+u)  +X8]  ~ 

^ooi'®'  - W ^ooo'*' 


where  h(8)  = (s+u)^  + 3>8  (s+u)  + 2>^s 


^000^*^  ■ 


8-{(8+U+>)  + 3X  (8+>)  (8+U+>)  + 2>  (s+X) 


- 30  - 


This  does  not  seem  to  reduce  to  simple  closed  form  expressions 
though  of  course  the  general  nature  of  the  solution  is  clear  in  terms  of 
the  roots  of  the  denominator. 

One  can  consider  the  limiting  form  (steady  state)  for  the  probabilities 
since: 

lim  s P«««(s)  = the  steady  state  value  , with  similar  expressions 

s-M  000  000 

holding  for  the  other  probabilities. 

Now  P ^(s)  = — + terms  corresponding  to  roots  of  the  cubic  expression 

uuu 

in  the  denominator  of  ^qqq{s).  By  setting  saO  in  the  algebraic  formulation  for 
the  partial  fraction  expansion  of  PqOO^*^  obtains 

a= ; r » where  p=  — 

1+30+6  p^+6 

Therefore.  p^Q^  = Urn  s Pqqq(s)  = a (as  given  above) 


and  similarly  Pqqj  = P PqoO' 

P012  ~ POOO’ 

Pl23  ~ POOO* 

These  same  values  will  be  obtained  as  a special  case  of  the  steady  state 
condition  in  Section  3.  1.  2.  2. 

It  is  interesting  to  note  that 

Pooo  + 5 Pool  + P0I2  + ^ Pl23  = = 1 ol  cour.., 

should  be  the  case. 

3.  1.  2 The  Steady  State 

The  general  steady  state  equations  are  obtained  by  setting  the  derivatives 
equal  to  zero  in  the  equations  of  Section  3.  1.  1.  1.  Steady  state  values  are  in- 
dicated by  a probability  notation  without  dependence  on  t,  e.  g. , Pqqq 
probability  that  no  components  are  in  the  service  system. 
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3»  1.2.  I The  General  Case 

When  each  of  the  three  components  are  allowed  to  have  their  own  values 
of  > and  U,  the  one  repairman  steady  state  equations  are  as  follows; 

(>1+>2+>3)  Pqoo  = ^1  Pool  ■*’  ^2  P002  * ^3  P003 

(Ui+>2+>3)  Pool  = POOO  *^2  P012  ^ ^3  ^013 

(u2+>i+^3)  Poo2  " ^2  Pqoo  ^i  P021  ^3  P023 

(U3+>  j+>2)  P003  “ ^3  Pqoo  ^1  Po3i  ^2  ^032 

{Uj+^3)  P021  " ^2  Pool  * ^3  P2I3 
(Uj+>2)  Po31  " ^3  Pool  ^ ^2  ^312 
(a2+>3)  P012  = ^1  P002  *^3  P123 

(U2+^l)  Po32  ~ ^3  P002  ^1  P32I 

(U3+>2)  P013  ~ ^1  P003  ^2  Pl32 

(U3+Xj)  Pq23  = ^2  P003  ^1  P23I 

^3  Pl23  " ^ 1 P023 
^3  P213  ' ^2  P013 
*^2  Pl32  " ^1  P032 

^2  P312  " S P0I2 
‘^l  P23I  ^2  P031 
^1  P32I  “ S P02I 

The  last  six  equations  yield  six  of  the  probability  values  directly  in 
terms  of  another  six  values.  One  can  write  the  other  ten  equations  in  terms 
of  ten  probability  values  by  using  the  relations  given  in  the  last  six  equations. 
This  procedure  results  in: 

()i+X2+>3)  Pqoo  “ ^1  Pool  ^2  P002  ^3  P003 

_ ^1  ^^2  ^3 

Pooi'ttj  P000'*’aj  Poi2‘''a]J'  Poi3 


- 32  - 


^2  ^3 

Po02  " oj"  PqOO  Po2l‘*'a^  P023 

^3  *^1  ^2 

P003  " ITj  PoOO‘‘’aJ  Po31  ■*■  Po32 


^Ui  +X3 

V >2  • ^021  " Pool  ^ Poi3 

+^2  ^ 

\“T^  J Po3i  " Pool  ^ P012 
/'^2  ■’’S  ^ _ 

\^)r^  / P012  " P002  P023 

/'  ^2  1 N 

\ r~  / Po32  “ P002  P021 


U-  + >. 


li-  + > 


■ ) P013  ■ P003  ■'■  P032 
" ) P023  * P003  ■*■  P031 


Where  = Uj  + X2  + >3,  a2  = U2  + Xj  + X^,  and  + >2  • 

In  principle  this  system  can  be  solved  in  algebraic  terms  to  yield  the 
probabilities  as  functions  of  the  six  parameters  u.,  >.,  i=l,2,3.  However, 
the  results  become  extremely  involved  and  it  is  not  considered  to  be 
practical  to  proceed  in  this  ' direct"  way. 

Numerical  solutions  of  the  linear  algebraic  equations  seem  to  be  the 
most  reasonable  way  to  obtain  particular  results.  Some  examples  are 
given  in  the  Appendix.  In  addition,  the  special  case  of  equal  components 
is  considered  in  the  next  section.  An  alternative  approach  to  utilization 
of  these  equations  for  the  study  of  component  availability  is  discussed  in 
Chapter  5. 

3,  1.2.2  All  Components  The  Same 

When  the  three  components  are  the  same  u.  = u and  X^  = X for  i=I,  2,  3. 
Moreover,  many  of  the  probabilities  are  necessarily  equal,  e.g.,  Pjqi  = 


L 
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P002"P003'  P012”P021"P013"P031  ®°  forth.  Typical  probabilities 

used  are  p^^^  (unique),  p^^^  (three  like  this),  p^^^  (six  like  this),  and 
P123  (six  like  this).  The  equations  defining  these  probabilities  are; 

' Pqoo"^  Pool 

(u+2»Poqj=>  Pqqq  + 2u  Pqj2 

(U  + ’)Poi2”’  Pool  ^Pl23 

^Pl23=  ^P012 

Let  P = ^ » then  the  solution  is; 

Pool  " ^ Pooo 
P012  “ ^^POOO 


T23 


= P P 


000 


Since  the  sum  of  all  probabilities  is  unity  one  finds; 

1 

Pooo  ^ j.  2 . . 3 

l + 3p+6p  +00 

These  results  agree  with  the  special  limit  case  given  in  Section  3.  1.  1.  2. 
3.  2 Two  Repairmen 

When  there  are  three  components  and  two  repairmen,  one  must  consider 

a total  of  ten  service  system  probabilities  P...  (t)  where  k and  j specify 

ijk 

components  in  service  and  i specifies  a component  waiting  for  service  with 
both  repairmen  occupied.  The  quantities  to  be  considered  are  (all  are 
fractions  of  time  in  the  general  case):  p^^Q,  p^^j,  p^^^.  Poo3'  Poi2'  Poi3' 
Po23'  Pl23'  P213'  P3I2'  indices  specifying  service  are 

unordered,  they  simply  tell  which  one  or  which  two  components  are  in 
service. 

3.  2.  1 The  differential  equations 

The  usual  procedure  is  followed  in  this  section.  Values  of  the  service 
probabilities  at  a time  t+At  are  expressed  in  terms  of  appropriate  values  at 
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time  t.  Higher  order  terms  in  bt  may  be  omitted  in  the  formulation,  or  if 
included  as  part  of  the  initial  concept  stage,  fall  out  upon  forming  limits 
as  At-0.  Somewhat  more  care  is  required  in  the  formulations  for  this  case 
than  in  the  cases  delt  with  in  previous  sections,  as  discussed  below. 


3.  2.  1.  1 The  General  Case 

Because  of  the  two  servicemen,  it  is  possible  (as  a lorical  case)  that 
two  distinct  service  activities  could  be  completed  in  the  same  time  interval 
^t.  Though  a logical  possibility,  this  is  in  fact  an  event  of  higher  order 
when  using  independent,  exponential  service  time  distributions.  Due  to 
the  somewhat  complicated  nature  of  this  case,  an  initial  formulation  of 
the  equations  is  given  below,  followed  by  the  resulting  system  of  differential 
equations.  Some  higher  order  terms  are  included  within  the  initial  form- 
ulation. 


P„„(,(t+»|=[l-()i+)2t>3)i»lPooo<*>^“l“PoOlW+“2‘'Po02l*>+“3"^003“' 

+ higher  order  terms  in  bt. 

PpOi(t+frt)=(l->  , 

a2^t(l-u^At)(l-X3it)  ^ ^3**  ^013^^^ 

U3At  Po23^*^ 

PQQ3(t+At)=(l-XiAt)(l->2At)(l-U3At)P^,03(t)+>3At  PQoo(t)+UjAt  Po23(t) 

Pqi 2^^'*'  "^3  1 "Uj  if)(  1 ^^^^012^^^^^  2^^^^  1 ^001 

>jAt(l-U2At)  PQQ2(t)  + U3At  Pi23(t)  + U3At  P2,3(t) 

Pq  1 3 (t+  At)-^(  1 - > 2 ( 1 -U 1 At) ( 1 -U3  At)PQ  j 3 (t)+>  3 At(  1 -u  j At)  Pqq  j (t)4 

>jAt(l-U3At)  ^003^^^  ^2^^  ^312^*^  ^2^*  ^123^*^ 

P023(t+it)=(l->iAt)(l-U2At){l-U3At)PQ23(t)  + >3At  (l-U2At)  Poo2<*>^ 

> At  <^■^3^^)  ^003^*^  *^1^*  ^213^*^  + UjAt  P3i2^^^ 

Cl 

Pj23(t+ At)=(l.U2At)(l.U3At)  Pj23(t)  + > ^ At  Po23(t) 

P2j3(t+At)=(l-UjAt)(l-U3At)  P2i3(t)  + Poi3(t) 

P312(t+At)=(l-UjAt)(l-U2At)P3j2  ^012 


Upon  taking  the  limits  as  At-0  after  first  dividing  by  At  and  forming 
appropriate  differential  quotients  one  obtains  the  following  system  of 
differential  equations  for  the  service  system  probabilities. 


This  system  is  complicated  for  a general  solution  in  terms  of  closed 
form  expressions  depending  on  the  six  parameter  values.  Special  cases 
are  considered  and  some  numerical  solutions  are  given  in  the  Appendix 


In  this  case  u.=u  and  >.  = > for  i=l,  2,  3.  It  is  also  true  that  Po01^*'^~^002^^^ 

correspond  to  one  in  the  service  system,  two  in  that  system,  and  three  in 
that  system  (with  one  waiting)  respectively.  Typical  quantities  are  used 
for  each  case,  the  probabilities  used  are  (t),  P»,  ^(t),  and  P.-,(t). 
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The  differential  equations  for  this  case  become; 


dt 

dt 

3t 

dP^23(t) 


= - Pooo'*'  + ^001 1*' 

= - Po01<*'  + ’ '’ooo'*'  + ^012''* 

= - (2u+l)  PoijIO  + 2'  Pool**'  ■''  ^ ^123*'* 

= -2uPj3,(t)t>P„j2(t) 


As  observed  i';  Sectic  a 3.  1.  1.  2 for  the  one  repairman  case,  it  does 
not  seem  possible  to  get  useful  direct  solutions  to  this  system.  Numerical 
values  can  be  computed  and  the  steady  state  limit  considered  (as  was 
done,  using  the  transformed  system  for  the  one  repairman  case).  It  is 
felt  that  these  techniques  are  sufficiently  well  illustrated  for  other  cases 
and  no  further  consideration  is  given  to  this  case  in  this  report. 


3.  2.  2 Steady  State 

The  steady  state  equations  for  the  three  component,  two  repairmen 
case  are  obtained  in  the  usual  way  from  the  differential  equations  of 
Section  3.  2.  1.  1.  The  notation  p..,  is  the  same  as  in  the  time  dependent 
case  except  that  the  service  system  probabilities  do  not  depend  on  time 
in  the  steady  state. 


3.  2.  ?.  1 The  General  Case 

The  linear  alpebraic  system  which  defines  the  service  system  steady 
state  probabilities  is  as  follows: 

(1  i + ^2'*'^3^  PqOO  " ^1  Pool  ^ ^2  P002  ^3^003 

(ai  + >2+>3)  Pqq2  = '^2  Pqqq  + U2  P012  ^3P013 

P002  " ^2  POOO  ^1  P012  ^ ^3  P023 

{U3+’j  + >2)  Poo3  ""  N PoOO  P013  ^2  P023 

(Ui-*-k'.2'’'’3^  P012  " ^2  Pool  P002  ^3  ^Pl23  P213^ 

(Uj+U3+>2)  Poi3  “ S Pool  P003  ^ ^2  ^P312  Pl23^ 

(u2+U3+‘‘  j)  Pq23  = ' 3 P002  ^2  P003  *^1  ^P213  P312^ 


{U2+U3)  Pi23  ~ ^1  ^023 
(U1+U3)  P213  = ^2  P0I3 
(Ul+i^)  P312  = ^3  P012 

It  does  not  seem  reasonable  to  use  this  system  to  express  tne  service 
system  probabilities  directly  as  functions  of  the  six  parameters  li.,  >. 
i=l,2,3.  Numerical  solution  of  the  system  for  various  parameter  values 
is  a much  more  reasonable  approach  to  the  generation  of  specific  results. 
Some  numerical  results  are  given  in  the  Appendix.  An  alternative  approach 
is  discussed  in  Chapter  5.  The  special  case  with  all  components  equal  is 
treated  in  the  next  section. 

3.  2.  2.  2 All  Components  The  Same 

In  this  case  UpU  and  >.  = > for  i=l,  2,  3 and  p^^ j=Pqq2“Po03’  f^012“f^013 

=Po23.  and  Pi23=P2 13^^31 2’  probabilities  represent  one  in  the  service 

system,  two  in  service  system  and  three  in  service  system  with  one  waiting 
for  service.  The  quantities  Pqqj,  PqiZ’  ^123  used  as  typical 

quantities.  The  governing  equations  for  this  case  are: 

^ Pqoo  ~ ^ Pool 

(u+21)  Pqqj  = > Pqqq  + 2u  P012 

{2a+>)  Pq^2  = Pqqj  + 2u  Pj23 


P012 


The  solution  is: 


Pool  " ® Pooo 
P012  “ Pooo 
Pl23  = ^ P ^ PqOO 

where  there  are  three  probabilities  of  each  type  and 


Pooo 


1 + 30+3  0^-+  -^p- 


I 
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Some  general  results  in  the  case  of  all  components  equal  are  given  in 
Chapter  4. 


3,  3 Summary  Results 

Because  the  solutions  for  general  parameter  values  are  so  difficult 
to  express  in  terms  of  and  summary  results  are  discussed  for  the 
all  components  equal  case  as  was  done  for  two  components  in  Section  2.  2. 


Three  components,  with  all  components  the  same,  allow  four  cases 
without  generalization  of  the  basic  queue  assumptions  (e.g.,  requiring 
FIFO  conditions).  These  are  denoted  by  Q®  q,  j,  2»  ^3  3* 

The  extreme  cases  ^ and  ^ treated  in  general  in 

Section  1.  2. 

Unlike  the  situation  for  two  components  in  Section  2.  2,  the  time  varying 
expressions  for  service  system  probabilities  and  hence  component  avail- 
ability have  not  been  obtained  for  the  cases  j and  2*  Some  numerical 

results  can  be  found  in  the  Appendix.  However,  the  steady  state  values  can 
be  summarized  here. 

S o 

Q steady  state  value  y =0 

^3,  1 " PqOO  ^ Pool  ^ Poi3 

(based  on  component  2 for  example) 


2 

1 + 2 0 4 2 p 

7~T  , 3 

l + 3p+6c  + 6p 


y 


“ Pooo  ^ Pool  P013 


Note  that  the  notation  differs  between  cases  and  2» 

latter  when  two  are  in  the  service  system  only  one  case  results  (both  in 
service)  in  the  former  two  cases  (since  either  customer  can  be  the  one  in 
service). 

2 

~2  1 + 2 p + p 

y z 3" 

l + 3p+3p  + 3p 

2 


31! 


yjjg^ 


1 

1 + P 


Direct  algebraic  comparison  establishes  that  ^ < y ^ < y ^ <y^for  all 
(positive)  values  of  p = — as  one  would  expect. 

17  2 3 

The  spread  between  values  like  ^ , y or  y , y serves  as  a potentially 
useful  measure  of  what  is  achieved  by  introducing  an  additional  serviceman 
into  the  service  system.  Figure  3 allows  such  a comparison  as  a function 
Qf  Of  course  one  repairman  gives  meaningful  values  of  availa- 

bility in  relation  to  the  zero  value  in  the  no  repairman  case.  An  additional 
repairman  is  seen  to  be  increasingly  worthwhile  as  p increases.  However, 
going  on  to  a third  repairman  produces  limited  effect  for  the  values  of  o 
shown.  As  p increases,  the  three  repairman  case  will  be  a real  improvement 
However,  such  values  yield  rather  low  availability  in  any  case  due  to  the 
tendency  of  components  to  fail  at  a faster  rate  than  they  are  repaired  (p>l). 
Such  component  availability  values  should  be  avoided  in  design  unless  there 
is  no  chance  to  avoid  them. 

Chapter  4 The  General  m Component,  r Repairmen  Case  when  All 
Components  are  the  Same. 

The  cases  of  two  components  with  one  repairman,  three  components 
with  one  repairman,  and  three  components  with  two  repairmen  have  been 
treated  in  previous  sections  for  all  components  the  same.  In  this  chapter, 
the  all  components  the  same  cases  are  given  general  treatment  for  m 
components  and  r < m repairmen.  The  cases  of  no  repairmen  and  m 
repairmen  were  treated  in  Section  1.  2 by  means  of  the  fundamental  relation, 
equation  (1).  However,  the  combinatorial  queue  model  can  be  used  to  obtain 
the  result  for  m repairmen  as  well.  An  example  of  such  a development  for 
m=3  is  given  by  the  equations. 


dPo(t) 


= - 3>  Pglt)  + 3U  Pj  (t) 


dPj(t) 

dt 


= - (U  + 2>)  Pj(t)  + > PQ(t)  + 2u  P2(t) 
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dP,(t) 

-gf—  = - (2u  + >)  P2(t)  + 2>  Pj(t)  + ^ P3(t) 


dP3(t) 


3u  P3(t)  + 3>  P^  (t) 


where  P.(t)  = probability  that  there  are  i components  in  the  3 repairman 
queue  at  time  t.  The  basic  assumptions  of  the  combinatorial  queue  model, 
previously  specified,  should  be  kept  in  mind. 

The  steady  state  values  are  determined  by  the  following  system  with 

r^'k  / li. 

Pr^Po 

(u+2  >)  pp)  Pq  + 2u  p^ 

(2u+>)  P2=2^  Pj  + It  P3 

P3=PP2 

2 3 1 

with  solution  p,  = r p^,  p2=c  Pq,  P3=p  p^  and  Pq= 2 ^ 


1 + 3 p + 3 p + p' 


The  availability  is  y ^ = pQ+2pj+p2  = ^ 3 “ 

1+3  pf  3 p + p^ 

in  agreement  with  the  result  from  Section  1.2. 

As  a general  approach  to  the  same  components  case  one  considers 

P|^(t)  as  the  probability  that  there  are  k components  in  the  service  system 

at  time  t,  out  of  a total  population  of  m components.  The  steady  state 

expression  is  p,  . In  the  case  of  Q ,p  indicates  k receiving  service 
K lYky  r K 

and  some  waiting  if  k<^r,  otherwise  r are  in  service  and  k-r  are  waiting. 
The  availability  for  a component  is  the  probability  that  the  component  is 
not  in  the  service  system.  If  k are  in  the  service  system  and  m-k  are  in 
the  remaining  (available)  population,  the  probability  that  a particular 
individual  component  is  available  is  (m-k)/m  on  the  assumption  that  all 
the  components  are  the  same.  Thus: 
y=Pr  (component  ' a''  is  available) 
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_ Pj.  (I  a'  is  in  non-service  population  | population  size  is  m-k) 

k=0 

Pr  (population  size  m-k) 

_ ^ (m-k)  ^ 

ks:0  m ^ 

Consider  the  case  Q®  , (m=2,  3 have  already  been  discussed),  for 

m,  1 

which  the  defining  equations  are: 

Pq  (t+&t)  = (l-m>M)  pQ(t)  + Uftt  Pj(t) 

It+ftt)  = [l-(m-k)  >&t]  (1-uAt)  P^(t)  + uAt  Pj^^j(t)+(m-k+l)  > At  P^.^U) 

for  k=l,  2,.  . . ,(m-l) 

Pn^(t+At)  = (1-uAt)  P^(t)  + >At  Pjn.i{t) 

These  yield  the  difference  differential  equations: 

dPo(t) 


dt 

dPk(t) 

3r 

dPit) 

m 

dt 


= - m>PQ(t)  + uPj(t) 


= -[(m-k)  >+u]  Pj^(t)  + uPj^^j(t)  + (m-k+l)>Pj^_j(t),  k=l,  2, . . . , (m- 1) 


= - u P^(t)  + > P^.i(t) 


The  steady  state  equations  are: 

Pi  = mp  Pq 

p^^j  = [l+(m-k)p]  pj^  - (m-k+1)  p Pj^_ ^ , k=l,  2,.  . . ,(m- 1) 

^Tn  ” ^ ^m-l 

Where  o=l/ii  as  usual. 

Induction  yields  the  solution: 
ml  .k  1-19 

Pk“  T^k)i  P Pq  ' •••» 


m 


and 


Po= 


1 


ml  k 
P 


m 
k=l 

This  result  can  be  used  to  obtain  the  m=2  and  m=3  results  with  one 


repairman  obtained  previously.  Thus 


I w 

\ > 


saiywatfrv 
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m=2  , Pq= 


1 + 2o  + 2p 


1 


”1=3  , p_  = n 5- 

l + 3e+6o  +6o 


The  availability  values  may  also  be  considered, 

m=2  . y=  r Pk  = Pq  ^ Pi 

k=0  ^ k u i 


y ~ y , as  previously  obtained  in  Section  2.  1.  2.  2, 

1 + 2 0 + 2 p 


m=3 


/=  r 


k=0 


^ 1^-  - Pk  = Po  + 4 Pi  3 P2 


y=  Pq  + 1 Po  3 Pq 


1 + 2 0 + 2 p 

y= Y 


, as  previously  obtained. 


l + 3p  + 6p  + 6p' 

Now  consider  the  cP  - case  which  only  has  meaning  for  m>  3 (except 
m,  c 

for  the  full  repairman  situation  when  m=2  which  is  not  under  discussion 
here).  The  service  situation  can  have  one  in  service  and  none  waiting, 
or  two  in  service  with  zero  or  more  waiting,  or  of  course  none  in  service. 
The  equations  governing  this  case  follow: 

Pj,(t+At)  = (l-m>  bX)  PQ(t)  + uM  Pj(t) 

Pj(t+ftt)  = [l.(m-l)  X /^t]  (l-uAt)  Pj(t)  + mXAt  P^CO  + 2uAt  P2(t) 


Pj^(t+it)  = [l-(m.k)XAt]  (l-2uAt)  Pj^{t)  + {m-k+l)X  At  Pj^_  j{t)+2uAt  Pj^^j 
for  2 < k < (m  - 1) 

P (t+At)=  (l-2uAt)  P_(t)  + XAt  P^  ,(t) 
m m m - 1 

This  leads  to  the  difference -differential  equations: 

dPo(t) 


(t) 


= - ml  Pn(t)  + a P,(t) 
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2 

I T 1 2 1+2  0 + P • 1 

= Pq  + 2p  p + p Pq  = ■ -j  — , as  previously 

1+3  p+3  p +^  p 

obtained. 

For  general  r the  situation  is  somewhat  more  involvec  tha.i  for  the 
r=l,  2 cases  shown  above.  The  major  feature  contributing  to  increased 
complexity  is  the  need  for  more  individual  equations.  The  defining 
equations  are: 

Pq  (t+&t)  = (l-m\At)  Pglt)  + uAt  Pj  (t) 

Pj  (t+At)  = [l-{m-l)XAt](l-uAt)Pj(t)  + mXAt  PQ(t)  + ZyAt 

P2  (t+At)  = [ l-(m.2)>At](1.2u  At)P2(t)  + (m-l)>At  Pj(t)  + 3uAt  P3(t) 

P^_  j(t+At)=[l  -(m-r+l)>  At]  [l-(r-l)  uAt]  Pj._i(t) 

+ (m-r+2)X  At  P -(t)  + ruAt  P (t) 
r - c.  r 

Pj^(t+At)  = [1 -(m-k))  At](l-ruAt)  Pj^(t)+(m-k+l)X  At  Pj^  j(t)+rUAt 

for  k=r,  r+1 , . . . , m - 1 

P^(t+At)  = (1-ruAt)  P^(t)  + >At  P (t) 
m m m- 1 

These  yield  the  difference-differential  equations: 

= - mX  PqU)  + u Pj(t) 

= -[(m-1))  + u]  Pj(t)  + mX  PQ(t)  + 2u  P2(t) 

= .f(m-2)>+2u]  P^(t)  + (m.l)XP^(t)  = 3u  P3  (t) 


dP^lt) 

~3r 

dPj(t) 
■ dt  ■ 

dP2(t) 

~^t“ 


p.ii..j'ii..  u 
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dP  i(t) 

■^—  = - r(m-r+l)>+(r-l)u]  Pr.i{t)  + {ni-r+2)X  Pr.2(t)  + ’’1^ 

dP,  (t) 

— ^ = - r(m-k)  X+ru]  P^(t)  + (in-k+l)>Pj^_j(t)+  ru  P,^^.i(t) 

k=r,  r+I , . . . , m - 1 


dP  (t) 
m' 

dt 


ru  p_(t)  + y P^_i(t) 


m 


The  steady  state  equations  are; 

Pi  = mP  Pq 

2p2=[(m-l)  0+1]  Pj  - mo  Pq 
3pj=[(m-2)  0+2]  p2  - (m-l)  p pj 


rp^=[(m-r+l)  0 + r-l]  p^_^  - (m-r+2)  p p^  ^ 

rp|^^j  = [(m-k)  o+r]  pj^  - (m-k+1)  OPi^_l  » k=r,  r+1 , . . . , m- 1 
rp  =0  p__  1 

This  system  does  not  lend  itself  to  a simple  general  solution  in  terms 
of  r.  One  can  develop  a solution  by  induction,  however,  it  does  not  seem 
to  be  worthwhile  in  view  of  the  resulting  complexity. 

Particular  cases  of  interest  can  of  course  be  treated  as  needed  from  the 
general  equations.  For  example  the  rs3  case  arrises  in  a four  component 
situation.  One  can  consider  the  r=3,  m=4  case,  for  which  the  equations 
become; 

Pi  = 4P  Pq 

2 p2  = (3o+l)  pj  - 4p  Pq 
3p3  = (2o+2)  P2  - 30  pj 
3P4=  PP3 


v*« ' v2>*u:. 
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^'hich  has  the  solution: 


Pi  = 4ppq 

P2  = 6pSo 

P3  = 4p^P0 

4 4 

P4=-J0  PO 


whe  re 


Po  = — 


l + 4p+6p  +4p  "^-^p 


The  availability  expression  in  this  case  is: 

Pk  = Po  Pi  +1  P2  +¥  P3 

,1  ,i2  , 3 l + 3p+3p^+p^ 

y = Po  + 3opo  + 3o  Po+P  Po  = 77rT3’:  4 4 

l + 4p4’6p  + 4p  ^ 

It  may  be  observed  that  y ^ which  is  the  full  repairman  case. 

This  should  be  the  case  since  there  are  only  three  repairmen. 

Chapter  5 Application  to  System  Design 


This  report  deals  in  detail  with  the  concept  of  component  availability 
within  the  context  of  various  system  service  procedures.  When  the  com- 
ponents operate  in  a system,  various  levels  of  availability  are  required 
for  individual  components.  It  is  therefore  of  value  to  consider  the 
different  component  availabilities  together  in  some  form  to  which  a trade 
off  analysis  can  be  applied.  Such  a formulation  is  presented  in  this 
chapter.  In  a particular  system  design  of  course,  the  required  level  of 
system  availability  drives  the  rest  of  the  availability  analysis.  System 
logic  is  utilized  to  relate  system  availability  to  component  availability 
levels.  The  consideration  of  this  chapter  relates  the  X and  U values 
(design  parameters)  and  component  availability  levels  for  two  component 
and  three  component  systems.  There  is  no  reason  that  the  procedure 
could  not  be  extended  to  more  components.  However,  the  underlying 
queue  problems  and  numerical  analyses  become  involved.  The  procedure 
will  be  introduced  for  a two  component  system  then  generalized  for  the 


■Via. 


! 
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three  component  case. 

5.  1 Two  Components 

In  a two  component  system,  two  component  availabilities  are  to  be 
considered,  they  are  denoted  by  y^  and  The  design  parameters  are 

the  failure  rates  Xj,  and  repair  rates  la^,  Because  of  the  difficulty 

in  combining  all  these  values  into  a relatively  simple  comparison  form  the 
following  procedure  has  been  developed. 

One  component  is  considered  to  be  the  reference  or  major  component, 
though  this  does  not  really  have  any  firm  significance,  and  the  other  com- 
ponent is  related  to  it  by  ratio  value.  In  addition  failure  is  related  to 
repair  by  a ratio  value,  assumed  to  be  the  same  for  both  components. 

These  assumptions  restrict  the  possible  cases  for  study  but  result  in  a 
very  managable  presentation  which  still  allows  a wide  range  of  options  . 

Following  the  procedure  outlined  above,  the  following  assumptions  are 

made,  taking  component  1 to  be  the  reference  (or  major)  component: 

u-.=0u,.  X,  = au,.  X-=au.,='^rPu, . One  can  consider  the  differential  equations 
2 '^’^1'  1 1’  2 2 1 

for  queue  probabilities  and  the  resulting  component  availabilities  as  functions 
of  time.  Only  detailed  computation  is  required  to  do  this.  However,  for 
purpose  of  illustration  only  the  steady  state  availabilities  are  treated  here. 

Making  the  substitutions  for  and  in  the  results  for  j steady 
state  (the  cases  q and  ^ great  interest  here),  one  obtains: 

^02  " ^ PqO 
Pol  ^ ® Poo 
Pl2  " ^ Poo 
P21  “ ° Poo 

i in  X M P X p\ 


where  A 3 aP 


For  fixed  0,  as  p-^  one  obtains  y 


IP=U)  = 2 1+3 

^ (a+l)^ 

5 for  P=0  which  agrees  with  the  numerical 


for  example  whena=l,  y^  — 
results  shown  in  the  figure 


l-«o  one  obtains 


which  shows  the  interchange  of  y,  and 


so  as  P-«>,  y^  and  y2- 

y^  under  the  transformation  P-  (by  comparison  with  the  P=0  case).  For 
example  as  for  a=l  one  finds  yj  = . 5 and  y^=  .25  (in  the  limit). 

In  this  formulation,  an  increase  In  0 means  that  failure  rate  is  increasing 
relative  to  the  ability  of  the  system  to  repair  components.  Of  course,  this 
reduces  component  availability.  For  a fixed  value  of  a the  availability  of 
component  1,  namely  y^  increases  with  increasing  P.  An  increase  in  P 
means  that  component  2 has  increasing  repair  rate  relative  to  component 
1 (when  P>1,  the  repair  rate  is  greater  for  component  2 than  for  component  1) 
Thus  it  is  not  the  repair  rate  that  produces  the  effect  on  availability  but  the 
waiting  time.  Even  though  the  repair  time  for  component  2 decreases  in  this 
instance,  its  availability  decreases  because  it  must  wait  while  component  i 
undergoes  service  with  relatively  longer  service  time. 

The  a,P  curves  of  Figure  4 can  be  used  in  a number  of  ways  in  design 
trade  off  studies.  For  example  if  system  availability  analysis  shows  that 
y > .75  and  y,>  .65  then  only  those  a,P  values  yielding  such  pairs  could 
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6.  1 Combinatorial  queues 


The  concept  of  a queuing  system  in  which  the  individual  members 
of  the  population  (components)  are  identified  and  followed  through  the 
stages  of  the  system  has  been  called  a combinatorial  queue  in  this  report. 

Such  queues  were  utilized  to  study  two  and  three  component  systems  with 
a FIFO  service  descipline  and  exponential  arrival  and  service  distributions. 

In  addition  to  more  general  service  times,  discussed  in  the  next  section,  a 
number  of  alternative  investigations  can  be  considered  for  the  study  of 
combinatorial  queues  as  models  for  availability  calculations. 

If  a component  can  be  replaced  with  some  short  (or  zero)  wait  then  one 
can  bring  in  the  concept  of  an  inventory  for  each  part.  The  down  part  may 
undergo  repair  then  enter  its  inventory  to  wait  until  needed  as  a replace- 
ment. This  sort  of  situation,  widely  studi  »d  in  the  literature  by  various 
methods,  has  not  been  treated  here.  However  it  can  be  incorporated  into 
a broader  description  of  the  combinatorial  queue  concept.  In  doing  so  it 
seems  likely  that  a more  general  definition  of  component  availability  would 
result.  One  way  to  do  this  would  be  to  define  the  component  as  the  set  of 
components  that  can  perform  a particular  role  within  a system.  Such  a 
general  component  is  otly  unavailable  when  the  set  becomes  empty  or  during 
the  transition  period  when  a failed  member  of  the  set  is  replaced  by  one  that 
is  operable.  There  seems  to  be  an  area  for  useful  research  in  this  direction. 

Another  potentially  useful  generalization  is  to  introduce  priorities  to 
replace  the  FIFO  type  of  service.  Different  components  can  be  assigned 
different  service  priorities  on  the  basis  of  their  importance  to  the  system  of 
which  the  components  are  a part.  These  priorities  would  represent  one  more 
system  input  to  the  component  availability  concept.  A closely  related  idea 
arises  when  there  are  several  repairmen  available  to  the  service  function. 
Repairmen  can  be  assigned  to  components  in  such  a way  that  when  a particular 
component  fails  it  must  be  repaired  by  a particular  repairman  (or  set  of 
repairmen).  An  extreme  case  of  this  for,  say  two  repairmen,  would  be  having 
one  repairman  dedicated  to  one  component  and  the  other  repairrpan  responsible 
for  the  repair  of  all  the  other  components.  Such  situations  are  likely  to  have 
useful  application  such  as  when  one  component  is  particularly  complex  or  when 
its  service  requires  special  safety  features. 
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Rather  than  have  each  component  distinct  in  its  descriptive  character- 
istics it  might  be  appropriate  to  have  several  components  the  same  forming 
a class  of  similar  components.  The  combinatorial  queue  would  then  deal 
with  a number  of  classes  where  the  particular  members  of  one  class  have  no 
special  features.  In  the  extreme  cases  where  all  components  are  the  same 
there  is  only  one  class.  This  case  has  been  treated  in  various  sections  of 
the  present  study.  In  many  respects  the  generalization  to  classes  of  compo- 
nents shows  potential  utility.  In  particular  it  might  lead  to  an  ability  to 
treat  larger  numbers  of  components  than  can  be  handled  directly  by  the  basic 
combinatorial  queue  models. 

6.2  General  Service  Times 


It  is  felt  that  a major  restriction  of  the  work  described  in  this  report  is 
the  use  of  exponential  ser\’ice  time  distributions.  A considerable  body  of 
research  has  been  carried  out  on  service  times  and  a number  of  distributions 
have  been  found  based  on  combinations  of  statistical  analysis  and  theoretical 
probability  arguments.  Therefore,  it  would  clearly  be  of  value  to  generalize 
the  concepts  of  this  report  to  a number  of  service  time  distributions.  The 
ideas  have  been  formulated  and  illustrated  here  for  the  exponential  case.  To 
generalize  one  would  have  to  employ  numerical  solution  techniques  and  simu- 
lation. It  is  unlikely  that  any  kind  of  "closed  form"  solution  approach  would 
be  possible. 

One  aspect  of  general  service  time  analyses  is  to  consider  the  service 
time  as  a sum  of  several  random  times  (possibly  dependent)  each  having  its 
own  characteristic  type  of  probability  distribution.  Those  characteristics 
arise  from  the  nature  of  the  service  being  carried  out,  possibly  depending 
upon  the  relative  amount  of  "search"  activity  involved  in  contrast  to  "repair" 
or  "replace"  activity.  Indeed  it  may  be  conjectured  that  service  time  distri- 
butions can  be  formulated  by  proper  mixes  of  such  basic  distributions.  Such 
studies  are  rather  distinct  from  the  ideas  of  this  report  but  combine  with 
them  in  a more  general  description  of  component  availability. 

Of  course  one  may  also  wish  to  generalize  the  way  in  which  components 
fail.  Exponential  failure  has  been  used  in  this  report.  This  kind  of  general- 
ization can  also  be  introduced.  However,  it  is  felt  that  the  service  time  is 
the  more  significant  effect  on  any  broader  study  of  component  availability. 
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6.  3 Waiting  Time  as  Input  to  Design  Analysis 

The  basic  resvilt  expressed  as  equation  (1)  in  Section  1.2  relates  compo- 
nent availability  and  waiting  time  in  the  repair  service  queue  for  components. 
Most  of  the  present  report  has  addressed  the  service  queue  directly  and 
attained  waiting  time  and  availability  for  special  cases.  An  alternative 
approach  is  to  consider  various  waiting  time  distributions  and  use  equation  (1) 
to  obtain  the  resulting  availability  values.  In  many  cases  this  wovild  be 
accomplished  by  numerical  integration  of  equation  (1 ) once  the  waiting  time 
distribution  was  specified.  Several  procedures  suggest  themselves  for  the 
formvilation  of  waiting  time  distributions: 


. Various  large  population  queue  problems  can  be  solved  to  give 
reasonable  limit  type  values  for  waiting  time  distribution, 

. Waiting  time  distributions  may  be  stipulated  in  either  niimerical  form 
or  as  specified  distribution  functions, 

. Results  for  simple  cases  with  known  solutions  for  waiting  time  can 
be  used  as  a basis  for  assuming  the  form  of  the  waiting  time  distribution  in 
more  general  cases. 


The  second  of  these  approaches  introduces  a synthesis  prdaleminto 
queuing  theory.  One  postulates  a desirable  waiting  time  distribution  and  uses 
equation  ( 1 ) to  obtain  the  corresponding  availability.  However,  it  is  necessary 
to  consider  whether  or  not  a service  queue  structure  exists  that  would  yield 
the  postulated  waiting  time.  If  such  a queue  or  queues  exist  the  next  problem 
is  to  construct 'an  appropriate  queue.  Theae  questions  of  queue  synthesis 
seem  to  hold  potential  for  considerable  amounts  of  useful  research. 


The  third  procedure  can  be  illustrated  by  considering  some  results 
obtained  in  this  report.  An  outline  of  the  argui'nents  involved  and  an  initial 
try  at  obtaining  the  form  of  a waiting  time  distribution  are  given  below. 
This  material  is  intended  only  to  illustrate  an  approach,  it  has  not  been 
developed  beyond  initial  stages. 


A large  class  of  waiting  time  distributions  may  be  considered  as 
steady  state,  not  depending  on  the  time  at  which  a component  enters  the 
service  queue.  For  the  two  components,  one  repairman  case  the  steady 
state  waiting  time  expression  was  found  to  be  of  the  form  ( for  FIFO  service): 


I 

a| 


1 
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w.  (0 ) = a.. 

W.  (s)  = 1 - b.  e’^^j  ® 

where  is  the  repair  rate  of  the  "other"  component  and  a™,  b,  are  nor- 
malized probabilities  of  none  in  the  queue  and  the  "other"  in  the  queue 
respectively.  Normalization  was  achieved  by  dividing  by  the  sum  of  the 
unnormalized  probabilities  which  represents  the  probability  of  component 
i not  in  the  queue,  that  is  the  availability  y.  of  component  i. 

One  may  generalize  from  the  above  result  and  consider  the  three  compo- 
nents, one  repairman  case.  In  that  case,  for  a typical  component,  say 
i = 1 one  may  consider; 


Vr  Pooo  + Poo2  + P003  + P023  P032  ' 


Wj  (0) 


PqOO 

^1 


where  A and  B depend  on  the  joint  distribution  of  service  for  components  2 
and  3. 

The  joint  density  for  two  independent  exponential  distributions  is: 


where  component  2 is  in  service  first  then  component  3. 

A reasonable  way  to  generalize  this  is  in  terms  of  the  moment 
generating  function  which  for  two  variables  is: 


I'.'i 


For  three  joint  exponential  variables,  supposedly  related  to  components 
1,  2 and  3 one  has  the  generating  function: 


This  yields  the  density  function: 


In  general  for  m joint  exponentials  the  density  function  has  the  form 


an  integro -differential  equation  for  component  availability  which  could  be 


BUujecteu  CO  analysis 


The  work  reported  on  in  this  report  has  delt  with  component  availability 
in  the  context  of  a repair  facility  present  as  part  of  the  total  operational  system 
The  point  of  view  considered  was  taken  from  the  operations  end  of  the 
system  and  the  desirability  of  available  components.  An  alternative  point 
of  view  is  in  terms  of  cost,  either  in  actual  exoenditure  or  modified  nerfnr. 


P ! 

I 


- 55 


mance,  of  desired  levels  of  availability.  The  effect  of  one,  two  or  more 
repairmen  has  been  investigated  as  regards  the  resulting  levels  of  component 
availability.  However,  the  cost  increments  have  not  been  considered.  For 
example  it  may  be  better  relative  to  some  measure  of  cost  (or  effective- 
ness) to  have  the  kind  of  availability  one  repairman  can  supply  rather  than 
have  two  repairmen  at  increased  expense.  One  can  formulate  an  appropriate 
objective  function  or  goal  of  system  operation.  These  qviantitative 
expressions  combine  system  requirements  in  terms  of  component  avail- 
ability and  the  cost  of  repair  operations.  They  should  also  include  the 
cost  of  availability  levels  representing  the  value  of  increased  availability 
and  the  penalty  for  reduced  availability.  A study  of  goal  seeking  or  opti- 
mlzStiXOZl  of  offocfXVw  CwJsCwXVOS  C&H  bs  to  tllO  COHStSTMlxivS 

imposed  by  relations  between  component  availability  and  the  repair  function. 
Such  relations  have  been  formulated  in  the  present  work  so  that  the  initial 
steps  have  already  been  taken  toward  trade  off  analyses  of  available  systems 
in  terms  of  component  availability.  Every  system  has  its  own  particular 
characteristic  so  that  a general  theory  of  trade  off  optimization  can  not  be 
made.  However,  some  general  methodological  procedure  can  be  formulated 
and  illustrated  using  special  system  configurations. 
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Appendix 


Some  typical  numerical  results  are  presented  here  in  graph  form. 
These  are  meait  to  be  illustrative  of  the  kind  of  results  one  may  . 
obtain  from  both  the  differential  equation  systems  and  the  steady  state 
formulations.  Typical  parameter  values  have  been  used  as  specified  on  the 
various  graphs. 

Differential  equation  systems  are  solved  using  a subroutine  recently 
developed  by  R.  Flynn  at  the  Polytechnic. 

The  systems  of  linear  algerbraic  equations  representing  steady  state 
values  are  of  course  redundant.  Thus  one  must  utilize  a subset  in  each 
case  and  our  procedure  was  to  omit  the  initial  equation  in  each  formulation. 
This  worked  very  well.  A standard  subroutine  from  the  Scientific  Sub- 
routine Package  was  employed.  Similar  calculations  were  used  to  produce 
Figures  4,  5,  and  6 of  Chapter  5. 

The  following  examples  are  presented  in  the  Appendix: 


Figure  Number. 


Case. 


2 time  variation 

Q . steady  state 
* 

Q . time  variation 
1 

Q . steady  state 

2 
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